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Abstract. As a component of the Flash Center's validation program, we compare 
FLASH simulation results with experimental results from Los Alamos National 
Laboratory. The flow of interest involves the lateral interaction between a planar 
' Ma=1.2 shock wave with a cylinder of gaseous sulfur hexafluoride (SFe) in air, 

and in particular the development of primary and secondary instabilities after the 
passage of the shock. While the overall evolution of the flow is comparable in the 
simulations and experiments, small-scale features are difficult to match. We focus 
on the sensitivity of numerical results to simulation parameters. 
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1. Introduction 



in 
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"^h \ The impulsive acceleration of a material interface can lead to complex 

Or 

^ fluid motions due to the Richtmyer-Meshkov (RM) instability. Here, 
H 

! the misalignment of pressure and density gradients deposits vorticity 



along the interface, which drives the flow and distorts the interface. At 



later times the flow may be receptive to secondary instabilities, most 



prominently the Kelvin-Helmholtz instability, which further increase 
the flow complexity and may trigger transition to turbulence. 

Verification and validation are critical in the development of any 
simulation code, without which one can have little confidence that 
the code's results are meaningful. FLASH is a multi-species, multi- 
dimensional, parallel, adaptive-mesh-refinement, fluid dynamics code 
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for applications in astrophysics p. Calder et al. discuss initial valida- 
tion tests of the FLASH code [2j • Herein we continue our validation ef- 
fort by comparing FLASH simulations to an RM experiment performed 
at Los Alamos National Laboratory [Bill]. 

2. Experimental Facility and Data 

The experimental apparatus is a shock-tube with a 7.5 cm square cross- 
section, as shown in Fig. ^ Gaseous SFg flows from an 8 mm diameter 
nozzle in the top wall of the shock-tube, forming a cylinder of dense 
gas in the otherwise air-filled test section. A Ma=1.2 shockwave travels 
through the shock-tube and passes through the cylinder. Our interest 
is in the resulting evolution of the SFg. All the experimental data is 
obtained in a plane normal to the cylinder axis, 2 cm below the top 
wall of the test section. The experiment is nominally two-dimensional; 
however, air diffuses into the SFg column as it flows downward, thick- 
ening the interface and reducing the peak concentration of the heavy 
gas. 

The initial SFg distribution (before the shock impact) is visualized 
directly by Rayleigh-scattering from the SFg molecules [3]. The pixel 
intensity in the experimental image gives only the mole fraction of 
SF6 relative to the peak mole fraction, Xsf6, which must be assumed. 
The distribution of SFq is only approximately radially symmetric, and 
the signal is dominated by noise at the level of about 5-10%. Smooth 
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Test Section 




Figure 1. LANL shock tube dimensions with test section detail. 



initial conditions for our simulations are obtained by fitting a radially- 
symmetric function to the experimental data. 

During the experiment the SFg distribution is indirectly visualized 
by visible light scattering off water /glycol droplets, which are seeded 
in the SFg. A sequence of experimental images is shown in Fig. [21 The 
shock traverses the cylinder in less than 25 /xs. The vortex Reynolds 
number of the flow, as measured experimentally, is Re= T/v ~ 5 x 10 4 , 
where T is the circulation and v is the kinematic viscosity. Each image 
is taken from a different experimental run. The water /glycol droplets 
can also be used to construct two-dimensional velocity vectors in the 
image plane using particle image velocimetry (PIV) 0, but the entire 
test section must be seeded, so simultaneous velocity and composition 
measurements cannot be obtained. 



hedla04_proc.tex; 2/02/2008; 11:04; p. 3 



4 



Weirs, Dwarkadas, Plewa, Tomkins, Marr-Lyon 



CC( CX.: 

Figure 2. Experimental time series of the distribution of SFe- The first image corre- 
sponds to 50 fis after shock impact; following images are at 190, 330, 470, 610, and 
750 fis. Intensity corresponds to the mole fraction of SF6. 

3. Flowfield Evolution 

As the shock traverses the cylinder, vorticity is deposited along the 
interface due to the misalignment of the pressure gradient (normal to 
the shock) and the density gradient (normal to the interface.) The 
density gradient arises from the gas composition; SFg is about five 
times as dense as air. Once the shock has passed through the SFg, 
the flow is dominated by a counter-rotating vortex pair, as shown in 
Fig. |2 Instabilities develop along the distorted interface at the edge 
of the primary vortices. The development and evolution of the vortex 
pair and subsequent instabilities at the interface proceed in a weakly 
compressible regime. More precise descriptions can be found in the 
references 0011]. 

The flowfield evolution is driven by flow instabilities and vortex 
dynamics, which are sensitive to the initial conditions and noise in 
the system. For validation this sensitivity is desireable because it pro- 
vides a severe test for the FLASH code. Figure El shows a sequence of 
images from our baseline simulation. The minimum grid resolution is 
78 microns, the initial peak mole fraction of SFq is 0.6, and the Courant 
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50 190 330 470 610 750 

Figure 3. Evolution of the SFg, with time elapsed after shock impact listed in fis. The 
mass fraction of SF6 is shown from a simulation in which Xsf6=0.6 and CFL=0.8. 

(CFL) number is 0.8. Overall the flow features in the simulation results 
are similar to those in the experimentally obtained images. Next we 
describe the effects of several simulation parameters on the computed 
results. The amount and location of small-scale structure, relative to 
the experimental data at 750 /j,s, will be used as a qualitative metric. 



4. Results 

In our initial investigations we have focused on the sensitivity of the 
computed solutions to several simulation parameters. We have consid- 
ered the dependence on the initial maximum mole fraction of SFg, the 
mesh resolution, the mesh refinement pattern, and the Courant number. 
We have also compared velocity data, and are beginning to consider 
three-dimensional effects. Here we show results only for different mesh 
refinement patterns and Courant numbers. More thorough discussion 
of the results can be found in [HJ. 

Regarding the initial maximum mole fraction of SFg, we find that 
simulations where Xsf6 = ®-6 seem to match the experimental results 
better than when Xs^6=0.8. At the higher value, the initial density gra- 
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dient is larger; this leads to greater vorticity deposition, faster instabil- 
ity growth, and consequently, excessive small scale structure. However, 
the time sequences match better for Xsf6=0-8- 

It is known that unavoidable discretization errors at discontinuous 
jumps in grid resolution can act as sources of spurious small-scale 
structure [Sj- To test this possibility we have run simulations in which 
a predetermined area around the vortices is uniformly refined to the 
highest resolution. Compared to fully adaptive refinement (the default) 
this approach significantly reduces the amount of perturbations intro- 
duced by the grid adaption process but increases the computational 
cost of the simulations. 

In Fig.|3]we compare the results from a fully adaptive grid and grids 
with maximally refined rectangles of 3 x 3 cm, 4x4 cm, and 4x8 cm. 
The vortex structure is always less than 2 cm across. For the different 
grids the large scale morphology remains the same, but the shape of the 
cross-section visibly differs depending on the grid, as does the amount 
and location of small-scale structures. In particular, differences are 
noticeable in the small-scale instabilities present on the vortex rolls. 
Since all other simulation aspects are the same, the differences must 
originate with perturbations at jumps in refinement. 

We then repeated the simulations on the different grids, but at a 
limiting Courant number of CFL=0.2. The results are shown in Fig- EI 
We observe much less variation between solutions on adaptive and lo- 
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Adaptive 3x3 rect 4x4 rect 4x8 rect 

Figure 4- Solutions on different grids, 750 /is after shock impact at CFL=0.8. Left 
to right: fully adaptive grid; 3x3 cm refined rectangle; 4x4 cm refined rectangle; 
4x8 cm refined rectangle. In the rightmost image, the refined rectangle covers the 
entire spanwise extent of the test section. 




Adaptive 3x3 rect 4x4 rect 4x8 rect 

Figure 5. Solutions on different grids for CFL=0.2; for other details see Fig. |1] 

cally uniform grids at CFL=0.2 than at CFL=0.8. One explanation for 
these results is that the errors at the fine-coarse boundaries are larger 
and lead to stronger perturbations at higher CFL numbers. An alterna- 
tive explanation might be that at higher Courant numbers, PPM does 
not adequately compute solutions at these conditions. Our simulations 
indicate that for FLASH, a lower CFL number leads to more consistent 
results on different grids. 

In addition to the SF6 distribution, experimental measurements of 
the velocity field in the vicinity of the vortex pair are available. The 
particle image velocimetry technique provides two-dimensional velocity 
vectors in the image plane. We find that in the frame of reference of 
the vortices, the velocity magnitude and large scale structure of our 
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simulations match those measured. The greatest discrepancies are at 
smaller scales at the outside edges of the vortices. 

Finally, we executed a speculative three-dimensional simulation. 
Our extension of the initial conditions in the third (cylinder axis) 
dimension is purely ad hoc, because we have no corresponding ex- 
perimental data. For this reason this simulation cannot be used as 
a validation test for the FLASH code, but we hope it will open a new 
line of investigation and discussion. In the cylinder-axis dimension, we 
varied the maximum mole fraction from X,s^6 = 0.64 at the top wall of 
the test section to Xsi?6=0.47 at the bottom, and in the same direction 
the radius of the cylinder increases slightly. These changes result in 
greater vorticity deposition at the top wall, and the vortex pair and 
instabilities evolve more quickly there. This behavior is expected based 
on the two-dimensional simulations; more interesting is that the maxi- 
mum flow velocity in the axial direction is greater than half that in the 
spanwise direction by the end of the simulation (750 [is). The stronger 
vortices have lower core pressures, and the pressure gradient in the 
vortex cores accelerates the flow from the bottom wall toward the top 
wall. The air is preferentially accelerated because of its lower molecular 
weight and confinement by the SFq, which acts like the wall of a tube. 
For the initial conditions we have assumed, the axial velocities suggest 
the three-dimensional effects are present. 
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To date we have made a large number of two-dimensional simulations 
to validate the FLASH code for problems dominated by vortex dynam- 
ics. So far, we have gained a better understanding of the sensitivity 
of the computed solutions to simulation parameters such as resolution, 
CFL number, and mesh adaption. While we can recover the overall 
morphology, the approximate amount and location of small scale struc- 
ture, and velocity field, we must make assumptions (though reasonable) 
about the initial conditions to do so. 

We continue to work on several fronts. We lack quantitative, phys- 
ically meaningful metrics for comparing simulation and experimental 
data. These metrics are difficult to develop and are rarely given the 
attention they deserve. FLASH simulations do not currently include a 
physical model for viscosity, but resolution-dependent numerical vis- 
cosity is present. Simulations with a minimum grid spacing of 78//m 
exhibit approximately the same amount of small scale structure as 
seen in the experimental data, while results on coarser grids show too 
little and on finer grids too much. We will soon begin simulations with 
models for physical diffusion; all the results presented here will then 
be reviewed. Our three-dimensional simulation, despite issues with the 
initial conditions, suggests that three-dimensional effects might be im- 
portant for this experiment. We are performing a systematic study of 



hedla04_proc.tex; 2/02/2008; 11:04; p. 9 



10 Weirs, Dwarkadas, Plewa, Tomkins, Marr-Lyon 

three-dimensional effects, and hope that experimental data will become 
available for comparison. 
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